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We analytically solve the Non-Markovian single electron spin dynamics due to hyperfine inter- 
action with surrounding nuclei in a quantum dot. We use the equation-of- motion method assisted 
with a large field expansion, and find that virtual nuclear spin flip-flops mediated by the electron 
contribute significantly to a complete decoherence of transverse electron spin correlation function. 
Our results show that a 90% nuclear polarization can enhance the electron spin T2 time by al- 
most two orders of magnitude. In the long time limit, the electron spin correlation function has a 
non-exponential 1/t 2 decay in the presence of both polarized and unpolarized nuclei. 

PACS numbers: 03.67.Lx, 72.25. Rb, 73.21.La, 85.35.Be 
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Spins in semiconductor nanostructures are promising 
qubit candidates for a solid state quantum computer 
because of their long decoherence times and potential 
scalability 0. To demonstrate the feasibility of a spin 
qubit, electron spin decoherence in semiconductor quan- 
tum dots (QD) has been widely studied both theoret- 
ically and experimentally 0, 0, 0]- The level 
discretization in a QD ensures that spin-orbit interac- 
tion induced spin relaxation is quite slow in QDs 0, 
which leaves the environmental nuclear spins, particu- 
larly abundant in III-V semiconductors such as GaAs 
(10 4 — 10 6 depending on the actual size of the QD), as 
the main source of decoherence for the electron spins. 

It has been shown that static thermal polarization of 
nuclear spins leads to inhomogencous broadening of elec- 
tron spins (which can be corrected usingthe spin echo 
technique 0) at a time scale of 10 ns and nuclear 

magnetic dipolar coupling leads to electron spin spectral 
diffusion and dephasing at a time scale of 10 ^s [9I Hc|. 
For the intervening period of time, the hyperfine interac- 
tion between the electron and nuclear spins can also lead 
to electron spin decoherence, which is in general non- 
Markovian because nuclear dynamics is slower than the 
hyperfine dynamics. 

The study of the non-Markovian electron spin dynam- 
ics in the presence of hyperfine interaction is a compli- 
cated problem due to its quantum many-body (1 elec- 
tron spin and N nuclear spins) nature, and has drawn 
wide spread attention recently 0, 0, 0, Q. Ana- 
lytically, an exact solution has been found in the case 
of a fully polarized nuclear reservoir |llj . while for the 
rest of the parameter regimes (in terms of nucle ar p o- 
larization and external field), perturbative theory [12j or 
effective Hamiltonians 0, ^3 have been used to study 
the problem. Numerically, only small systems with typ- 
ically less than 20 spins have been explored because of 
the extremely large Hilbert space [lj, UjJ Ha ■ 

In this Letter we focus on the problem of the spin deco- 
herence of a single electron due to hyperfine interaction 
with the surrounding nuclear spins. Although at a finite 



magnetic field the direct electron-nuclear spin flip-flop 
is highly unlikely due to the Zeeman energy mismatch, 
higher-order processes where electron spins do not flip 
are possible. For example, conduction electron medi- 
ated nuclear spin interaction (RKKY) has been stud- 
ied for a long time in both metals and semiconductors 
El • Here our focus is the backaction of the electron- 
mediated RKKY interaction between nuclear spins on 
the single mediating electron spin. We start from the ex- 
act electron-nuclear-spin Zeeman and hyperfine Hamil- 
tonian and use the equation-of-motion approach in the 
Heisenberg picture. Helped by a systematic large field 
expansion, we solve the full quantum mechanical prob- 
lem analytically and reveal the crucial importance of the 
electron-mediated nuclear spin flip-flop processes in the 
decoherence of an electron spin. 

Calculating Green's functions with the equations of 
motion is an old technique in solid state physics |l8j . 
What is novel in our current study is that we use this 
venerable technique to attack the new problem of spin 
decoherence, which is generally studied usin g q uantum 
master equations for the density operator |l9| . This 
traditional approach originating from quantum optics is 
more adapted in dealing with weak interactions between 
a system and its reservoir. Wc demonstrate in this study 
that a properly defined correlation function can be used 
to fully characterize the decoherence properties of a two- 
level system, and the equation of motion approach can 
be a powerful tool in studying non-Markovian dynamics. 

We model the coupled electron-nuclear-spin system by 
the Hamiltonian 



H = Lu S z + Y,A k I z k S* 
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where S and I represent electron and nuclear spin oper- 
ators respectively, loq is the external magnetic field, A k 
is the hyperfine coupling constant with the kth nucleus, 
and h = 1. In this Letter we assume I = \ for simplicity, 
though all calculations can be generalized for arbitrary 
/. For a two-dimensional QD with a Gaussian electron 
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wave function, A k has the simple form A k = Aoe~ k / N 
with k S (0, oo) [121 . For convenience we assume Aq = 1 
so that time is measured in the unit of l/A . 

To describe the decoherence between electron spin | f ) 
and | |) states, we introduce a retarded transverse spin 
correlation function 

G±(t) = -i6(t)(y o \S-(t)S + (0)\y o ) . (2) 

Here 9{t) is the usual step function, and is the 

initial wave function of the system where the electron 
and nuclear spins are assumed to be in a product state, 
with the electron having spin down initially, i.e. = 
| JJ,; I k , I k , • • • , I k ) |2£j . This spin correlation func- 
tion represents the phase fluctuations between electron 
spin up and down states in the presence of the nuclear 
spin reservoir, which can be most clearly seen in the 
Schrodinger picture 

G x (t) = -i9(m;I^,---,I§ N \e im / n S-e- iHt ^S + (0) 

= -i6(t){(it;I* ki ,---,I* kN \e im / h }s- 

x{e- im / n \t,I z kl ,---,I z kN )} ■ ( 3 ) 

The term in the first curly bracket represents the evo- 
lution of the electron spin down state in the presence 
of the hyperfine interaction, while the term in the second 
curly bracket represents the evolution of the electron spin 
up state in the same environment. If no electron spin 
flip occurs, any decay in the calculated average would 
be due solely to dephasing between the electron spin up 
and down states. Obviously, electron spin flip will also 
cause decay of the correlation function. Therefore, G± (t) 
contains the complete decoherence information for the 
electron spin in consideration. 

An iterative equation of motion (EOM) for the spin 
correlation function G± (t) can be obtained by differenti- 
ating G±(t) with respect to time and then perform the 
Fourier transform. In general, for two arbitrary operators 
A and B, 

uj((A: B)) u = <¥o|4(0)B(0)|*o) + (([A, H];B)) U , (4) 

where {(A;B)) U is the Fourier transform of 
(#oM*)-B(0)|*o>- For G±(t), we use G±(lj) to 
represent its Fourier transform, so that 

uG x (u>) = l + (([S-,H\;S + )) u . (5) 

The second terms on the right hand side of Eqs. and 
l|5|l involves the calculation of higher-order correlation 
functions. A cut off or decoupling scheme has to be ap- 
plied to eventually close the set of EOMs. 

After G±(lli) is obtained, real time dynamics of G±(t) 
can be easily calculated by an inverse Fourier transform 
using the spectral function defined as 

p(w) = -TmG±(w)/n. (6) 



In general there are two types of contribution to the 
spectral function after performing analytical continua- 
tion (lu — > lu + i0 + ) [J^: a delta function Z p S(uj — u} p ), 
and a non-vanishing imaginary part of the self energy 
resulting from branch cuts, which results from the inte- 
gration of continuous poles. The delta function leads to a 
coherent oscillation with a single frequency lu p , while the 
continuous part leads to dephasing in the time evolution 
of the spin correlation function G±(t). 

We consider the general case of partially polarized and 
unpolarizcd nuclear spin reservoir where both the num- 
bers of spin up and down nuclei are of order A. The 
difference in the numbers of the two spin species is char- 
acterized by an effective polarization P — (N^ — N±)/N, 
where (iVj ) are the number of nuclear spins in the up 
(down) states. Now the effective magnetic field takes the 
form fl = u>o + J2k A k (I k ), where (I k ) represents time- 
averaging of I k (t). We consider the physically relevant 
case of large effective fields (f2 ~ A, requiring that we 
have either a reasonably large external field, or a nuclear 
reservoir with finite polarization) , and focus on the spec- 
tral broadening near lu — £1, which leads to dephasing of 
transverse electron spin magnetization. 

Previous studies [111 |l2j indicate that when only the 
direct electron-nuclei spin flip-flop is considered, the de- 
cay amplitude of the electron spin correlation function 
is of the order 0(1 /N), and the correlation function has 
almost undamped oscillations. Clearly, such direct pro- 
cesses are energetically unfavorable in high effective mag- 
netic fields. However, if the higher-order virtual process 
(electron mediated nuclear spin flip-flop) is included, we 
expect that nuclear field fluctuation will give rise to com- 
plete decoherence in the electron spin. In other words, 
the delta function (indicating no damping) in the spec- 
tral function would be broadened (decoherence) after the 
virtual processes are included. The spectral weight in the 
low energy region where w ~ 0(1) has been found to be 
negligible H3|. 

In the following calculation we treat the nuclear field 
A k I k within the adiabatic approximation, which 
is physical since S~(t) has an oscillation frequency 
f2 ~ A while the nuclear field varies in a much 
longer time scale, so that in the Fourier transform 
of (*olE fc ^fcW<S" (t)5t(0)[* > we can simply re- 
place YlkAkI k {t) by J2kAk{Ik)- Gj_(uj) is related to 
the higher-order correlation function ((I k I k ,S~; S^)) u 
through the following equation 

- m) G ^ = 1+ i E A k A k >{(I k ll,S-;S^)) 

k=£k' 

(7) 

Here {(I k I k ,S~; S^))u represents nuclear spin k and k 
flip- flopping with each other while electron spin returning 
to its initial state (|). 

Calculating this higher-order correlation function re- 
quires a cut-off to terminate the iteration. The struc- 
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FIG. 1: Spectral function p(£o) (where u — lj — Q — N/8Q) 
with P = for Q = 2.5iV (left panel) and SI = 4.5JV (right 
panel). In both panels the solid lines represent the spectral 
functions calculated by considering only self-energy £i(<D). 
The dashed lines show the total contributions from Ei(w) 
and Y,2(6i)). The original delta functions S(u)) in the spectral 
densities are broadened due to the effects of the flip-flopping 
of nuclear spins. 



ture of the iterative equations reveals that a natural cut- 
off does exist. The correlation functions with one flip- 
flopped nuclear spin pair has two contributions to the 
sclf-cncrgy, one of them proportional to iV 2 /(4f2) 2 , the 
other N 3 /(Ail) 3 . For two pairs of flip-flopped nuclear 
spins the two contributions are proportional to -ZV 4 /(4S1) 4 
and 7V 5 /(4f2) 5 respectively. This geometrical series con- 
verges quite fast when Q > N (large field expansion). 
Physically, the expansion parameter N/Q appears be- 
cause the intermediate high energy state where the elec- 
tron spin is flipped requires energy SI, and N comes 
from the summation over all nuclear spins. In the limit 
of SI 3> N, only the first-order term of the self-energy 
contributes significantly. Neglecting iV 4 /Sl 4 and higher- 
order terms, the final expression of the spin correlation 
function is 



1 



(P 2 -1)7V 2 V /-v (1-P 2 )N3„ ,_V 
W TSnl 32Q3 2 -'2[^) 



(8) 
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with the lowest-order self-energy taking the form 



iCj 2 - 3)log 
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2a) - 1 



2a) + 1 



3|w| + l] , 



(9) 



for \Cj\ < 1/2. Here a) = u - O - N/8Q. The exact 
form of S 2 (o)) is also found [U- Both self-energy terms 
have branch cuts or non- vanishing imaginary parts when 
\uj\ < 1/2, leading to dephasing when calculating G±(t). 
Another significant feature of G±(u>) is that it does not 
have a (^-function component anymore, indicating that 
the decoherence of G± (t) will be complete. In addition, 
Eq. (jHJ indicates that the amplitude of G±(u>) is ~ O(l), 



FIG. 2: T2 determined from the half width of the peak of 
spectral functions for different nuclear polarizations (P) and 
effective fields (fi). In all cases the spectral functions are 
calculated using both £i(u>) and £2(01). We have chose N = 
10 5 &ndA = J2 A k= 92 ^eV. 



in contrast to the fully polarized case, where G±(cu) ~ 
0(1/N). 

Figure ^ shows the calculated electron spin spectral 
functions for different effective fields. We compare the 
results of including only Si (J)) (solid lines) and those 
with both Si (a)) and £2(0)) (dashed lines) for SI = 2.5iV 
and SI = 4.5./V. The two panels clearly show the validity 
of the large field expansion for SI > 2.5N. For smaller 
O more higher-order terms need to be included to at- 
tain convergence. Indeed, even if SI < N there is no 
divergence in our theory, since there could be at most 
Ni (Ni < iVj) flip-flopped nuclear pairs in the system, 
so that there is an upper limit to the number of EOMs 
and terms in self-energy. The right panel of Fig. ^ (f° r 
SI = 4.5./V) shows that the contribution of S 2 (cl>) is now 
completely negligible. Using hyperfine coupling constant 
of bulk GaAs [22j, we estimate that SI = 2.5A^ corre- 
sponds to a magnetic field of 5 Tesla. Figurc^also explic- 
itly shows that the original delta function in the spectral 
function is now broadened after taking into account the 
electron-mediated flip-flop of nuclear spins. According to 
Eq. ©, in the limit Q > JV or P = 1, both self-energy 
terms go to zero, so that the delta function form of 5(i2>) 
of the spectral functions would have been recovered, and 
there would have been no decoherence effect. 

The decoherence time T2 for the electron spin can be 
determined from the half-width (Aw) of the spectral peak 
(T-2 = 1/Ad)). Figure El shows T2 as functions of the 
nuclear spin polarization P and the effective magnetic 
field SI. It is clear that T2 only increases slowly with 
the external magnetic field, but is much more sensitive 
to the nuclear polarization. If the nuclear polarization 
P is raised to 0.9 from 0, T2 increases by almost two 
orders of magnitude. Physically this is quite reasonable, 
as increasing polarization would reduce the phase space 
for nuclear spin flip-flops, while increasing external field 
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FIG. 3: Time evolution of the envelope of Re{Gj_(i)} for 
various regimes of T2. G±(t) is obtained from the spectral 
function by inverse Fourier transform. It should be noted that 
the real evolution is modulated by the fast oscillation term 
e iSlt . We have assumed N = 10 5 and A = A k = 92 peV in 
these calculations. 



slowly reduces the cross-section of these processes. 

The real-time dynamics of G±(t) is obtained with the 
inverse Fourier transform Gj_(t) = —i6(t) J p(uj)e~' lu ' t duj 
using the spectral function calculated with both £i((D) 
and Y,2{u>). Figure|3plots the time evolution of the enve- 
lope of Re{G±(t)} for three different parameter regimes. 
The solid line represents the case of fast decay with 
il = 2.5N and no polarization. The dotted line shows 
that increasing the magnetic field can increase the co- 
herence time moderately. If the nuclei in the QD are 
polarized to 90%, the amplitude of the fast oscillation 
in electron spin (with frequency f2) could be maintained 
for a much longer time as indicated by the dashed line. 
Notice that here the amplitude of G±(t) does decrease 
as quickly as in the previous two cases initially, but elec- 
tron spin quantum coherence is only partially lost so that 
clear revival phenomenon is visible after even several /is 
in Fig. □ 

The long-time asymptotic behavior of G±(t) can also 
be extracted from p(u>), which is nonzero only when |t*>j < 
1/2. Calculating the inverse Fourier transform at t 3> 1, 
we find 



G±(Jk) oc 



144ft 2 



tt 2 (1-P 2 )N 2 t 2 ' 



(10) 



The 1/t 2 power-law decay here can be compared to the 
1/t power-law decay found in Ref. 0, for large mag- 
netic fields, where electron-mediated nuclear spin flip- 
flops are not taken into account, and the exponential de- 
cay found in Ref. 01 > where an effective Hamiltonian for 
the nuclear spin flip-flop is considered. 



In summary, we have presented a detailed analytical 
study of transverse electron spin decoherence using large 
field expansion. We find that electron-mediated nuclear 
spin flip-flops contribute significantly to electron spin 
dephasing by generating fluctuations in the Ovcrhauser 
field (the nuclear field) for the electron spin. We find 
that 80-90% nuclear polarization can enhance the elec- 
tron spin Ti time by two orders of magnitude into the fj,s 
time scale in a 5 T external field. We also show that the 
long time asymptotic behavior of the spin decoherence is 
1/t 2 . 
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